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ABSTRACT 



The intergalactic gas in groups and clusters of galaxies bears the indelible stamp of galaxy 
formation. We present a comparison between observations and simple theoretical models indicat- 
ing that radiative cooling governs the entropy scale that sets the core radius of the intracluster 
medium. Entropy measured at the radius 0.1r2oo scales as T^/"^, in accord with cooling-threshold 
models for the regulation of core entropy. Cooling of baryons to form galaxies is likely to lead to 
feedback, and the signature of feedback may appear farther out in the cluster. Entropy measured 
at the radius rsoo in all but the most massive clusters exceeds the amount that can be generated 
by hierarchical accretion. However, feedback that smoothes the density distribution of accreting 
baryons, perhaps via galactic winds, can boost entropy production at the accretion shock by a 
factor ~2-4. An initial comparison of entropy at r5oo to smooth accretion models shows that 
smooth accretion is a plausible explanation for this excess entropy and suggests that baryon 
accretion onto groups was smoother than baryon accretion onto clusters. 
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1. Introduction 

Clusters and groups of galaxies are interest- 
ing laboratories for studying the cooling and 
feedback processes that govern galaxy forma- 
tion because these non-gravitational processes 
have clearly altered the structure of intraclus- 
ter and intragroup media. Simulations show 
that the density profiles of dark-matter halos 
created by hierarchical structure formation are 
nearly self-similar, with a moderate trend for 
dark matter to be more concentrated toward 
the centers of lower- mass haloes (e.g., Navarro, 
Frenk, & White 1997). The baryonic density pro- 
files of those halos in purely gravitational sim- 
ulations are also nearly self-similar, leading to 
an X-ray luminosity-temperature relation that 
agrees with self-similar scaling: Lx cc T|- in the 
brcmsstrahlung-dominated regime (Muanwong et 
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al. 2001; Borgani et al. 2001). 

Cluster observations showing that Lx ~ 
for Tx ^ 2 keV were the first clear indication 
that galaxy formation somehow breaks the self- 
similarity of the baryonic component (e.g.. Edge 
& Stewart 1991; David et al. 1993; Markcvitch 
1998; Arnaud & Evrard 1999). This effect shows 
up even more dramatically in the steeper Lx- 
Tx relation for groups (Helsdon & Ponman 2000) 
and the surprising faintness of the unresolved 
~ IkeV background (Pen 1998; Wu, Fabian, & 
Nulsen 2001; Bryan & Voit 2001). Recent ob- 
servations resolving the temperature structure of 
clusters have revealed similarity breaking in the 
mass-temperature relation as well. Self-similar 
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models with oc T^' , where A/a is the to- 
tal mass within a sphere whose mean density is A 
times the critical density, agree well with purely 
gravitational simulations (e.g., Evrard, Metzler, 
& Navarro 1996; Bryan & Norman 1998). How- 
ever, X-ray observations that resolve the temper- 
ature structure of clusters now indicate a slightly 
steeper slope Ma ^ T^^ and a mass normaliza- 
tion ~1.4 times larger than in simulations without 
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galaxy formation (Horner, Mushotzky, & Scharf 
1999; Nevalainen, Markevitch, & Forman 2000; 
Finoguenov, Reiprich, & Bohringer 2001; Sander- 
son et al. 2003). 

Including galaxy formation in the simulations 
generally steepens the Lx-Tx relation, but it is 
not clear whether the dominant mechanism re- 
sponsible for similarity breaking is cooling or heat- 
ing. Simulations that preheat the intracluster 
medium by imposing an entropy floor equiva- 
lent to the 100-150 keVcm^ levels observed in the 
cores of groups by Ponman et al. (1999) produce 
good agreement with the observed Lx-Tx rela- 
tion (Bialek, Evrard, & Mohr 2001)."^ Yet, sim- 
ulations that include radiative cooling without a 
compensating heat source seem to agree equally 
well with the data, at the expense of an uncom- 
fortably large condensed baryon fraction (> 30%) 
on group scales (Muanwong ct al. 2001; Dave, 
Katz, & Weinberg 2002; Tornatore et al. 2003; 
Kay, Thomas, & Theuns 2003). 

The insensitivity of the Lx-Tx relation to the 
details of galaxy formation can be Tindcrstood in 
terms of the entropy threshold for cooling within a 
Hubble time (Voit & Bryan 2001). Gas at ~ 1 keV 
with entropy equivalent to < lOOkcVcm^ can ra- 
diate away its thermal energy in less than a Hub- 
ble time. It must either condense or be heated 
by feedback to a higher entropy level. The cool- 
ing threshold therefore imprints a physical entropy 
scale on a medium that would otherwise be nearly 
scale-free. A flattening of the entropy distribu- 
tion below the cooling threshold, not necessarily 
an absolute entropy floor, is all that is needed to 
produce the observed Lx-Tx relation, as long as it 
forces the cluster's density profile to be flatter than 
n(r) oc r^'^/^ where the cooling time is less than a 
Hubble time. Then, the bulk of the cluster's lumi- 
nosity comes from regions where the cooling time 
is similar to a Hubble time, linking the cluster's 
overall Imninosity to the cooling threshold (Voit 
et al. 2002). Simulations of cluster formation per- 
formed with a wide range of feedback efficiencies 
support this basic picture, as long as that feed- 
back heats the affected gas to > IkeV (Kay et 
al. 2003). In order to distinguish the effects of 

^Entropy in this paper is quantified in terms of the X-ray 

observable K = Tn^ ^ where is the electron density. 
The standard specific entropy scales with the logarithm of 
this quantity. 



cooling from those of heating, more information is 
needed, and some of that information can be found 
in observations that resolve the entropy profiles of 
clusters and groups (e.g., Ponman et al. 2003), 
which depend in detail on the interplay between 
cooling, heating, and hierarchical accretion (Voit 
et al. 2003). 

The purpose of this Letter is to highlight two 
trends found in those spatially resolved entropy 
observations and to present an initial compari- 
son with theoretical predictions by way of sim- 
ple analytical models. Section 2 focuses on the 
core entropy measured at 0.1r2oo, where r2oo is 
the radius at which A = 200. Early measure- 
ments of core entropy suggested that a global en- 
tropy floor ~ 135keVcm^ applied to all groups 
and clusters (Ponman et al. 1999; Lloyd-Davies, 
Ponman, & Cannon 2000), but updated measure- 
ments indicate that entropy at 0.1r2oo scales as 
r2/3 (Ponman et al. 2003), as expected if the 
cooling threshold indeed determines core entropy. 
Section 3 focuses on entropy measured at rsoo, 
where A = 500. Somewhat unexpectedly, entropy 
at these much larger radii also exceeds the predic- 
tions of self-similar models, sometimes by as much 
as ^ lO^keVcm^ (Finoguenov et al. 2002; Pon- 
man et al. 2003). Comparing these measurements 
to analytical models for smooth accretion suggests 
that smoothing of the intergalactic medium before 
it accretes onto a cluster or group is a promising 
mechanism for producing this entropy excess (Voit 
et al. 2003; Ponman et al. 2003). In § 4 we dis- 
cuss the implications of these results and offer a 
few hypotheses to test with numerical simulations. 

2. Core Entropy and Cooling 

The entropy of intracluster and intragroup me- 
dia is generated primarily by shocks owing to 
mergers and accretion. Simulations show that 
these gravitationally driven processes lead to an 
entropy profile K{r) ~ r^-^ (Borgani et al. 2001), 
a scaling that reflects the history of mass accretion 
onto the cluster (Tozzi & Norman 2001; Voit et al. 
2003). Alterations of that profile owing to non- 
gravitational heating and radiative cooling should 
be most obvious at small radii, where the gravi- 
tationally generated entropy is smallest. Models 
involving strong global preheating typically pre- 
dict that cluster cores should be nearly isentropic. 
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Fig. 1. — Relationship between core entropy and the 
cooling threshold. Each point with error bars shows the 
mean core entropy Ko,i, measured at 0.1r200! for eight clus- 
ters within a given bin of luminosity-weighted temperature 
^lum! a^iid small circles show measurements for individual 
clusters (Ponman et al. 2003). The dotted line shows a self- 
similar relation calibrated using the median value of i^o.i 
measured in simulation L50-(- of Bryan & Voit (2001) which 
does not include cooling or feedback. The solid line shows 
the cooling threshold Kc{T), defined to be the entropy at 
which the cooling time equals 14 Gyr, assuming the cool- 
ing function of Sutherland & Dopita (1993) for 0.3 solar 
metallicity. The dashed line shows the entropy at 0.1r2oo 
in the model of Voit & Bryan (2001) when this cooling func- 
tion is used. Solid points without error bars show the total 
radiative entropy loss -K^ad computed for individual ha- 
los with the accretion histories of Tozzi & Norman (2001, 
squares) and Voit et al. (2003, triangles). These simple 
treatments of cooling reproduce the slope but slightly un- 
derestimate the normalization of the observed -K^o.i(Tiuni) 
relation, which shows no sign of an absolute entropy floor. 



at the level of the global entropy floor (Balogh, 
Babul, & Fatten 1999; Tozzi & Norman 2001; 
Babul et al. 2002). In models that rely more 
on cooling to break self-similarity, the details of 
the core entropy profile depend on the interplay 
between cooling, feedback, and merging, and per- 
haps on conduction as well (Voit et al. 2002). Re- 
gardless of these details, the entropy scale of sim- 
ilarity breaking set by cooling should correspond 
to the cooling threshold KdT) oc T^/^A'^/^, which 
leads to Kc oc T^/^ when bremsstrahlung domi- 
nates the cooling function A. 

Recent observational evidence favors cooling- 



threshold models. Figure 1 shows measurements 
of core entropy JsTo.i drawn from a compilation of 
cluster observations. Ponman et al. (2003) col- 
lected data from 64 clusters and groups and pre- 
sented the mean entropy levels within temperature 
bins containing eight clusters each. Binning the 
data averages out the cluster-to-cluster variations, 
which can range up to a factor of three at a given 
temperature (e.g., Mushotzky et al. 2003). These 
measurements sit well above the dotted line show- 
ing a self-similar relation calibrated to match the 
purely gravitational simulations of Bryan & Voit 
(2001). However, the measurements do not level 
off at a well-defined entropy floor. Instead, they 
trace a power-law relation Ko,i oc T°-^^, nearly 
identical to the scaling of the cooling threshold 
(Ponman et al. 2003). 

Entropy profiles of individual clusters and 
groups likewise show no evidence for an abso- 
lute entropy floor. The cores of groups observed 
with XMM-Newton are not isentropic but instead 
have entropy profiles similar to those of clus- 
ters (Mushotzky et al. 2003; Pratt & Arnaud 
2003). In fact, scaling those profiles according 
to if oc T^/^ produces a much closer match to 
massive-cluster profiles than the expected self- 
similar scaling K oc T (Pratt & Arnaud 2003; 
Ponman et al. 2003). Furthermore, the core en- 
tropy levels required to explain the Lx-Tx rela- 
tion with isentropic models are far higher than 
those observed (Voit et al. 2003). 

Simplistic models involving the cooling thresh- 
old come close to matching the /^o.i-Tium relation, 
but the match is not exact. The solid line in Fig- 
ure 1 shows the cooling threshold Kc{T\u^), de- 
fined to be the entropy level at which the cooling 
time equals 14 Gyr, according to the cooling ftmc- 
tion for 0.3 solar metallicity from Sutherland & 
Dopita (1993). This threshold lies ~25% below 
the binned data points with a very similar slope. 
The dashed line shows -ftTo.i computed from Kc 
using the prescription of Voit & Bryan (2001; see 
also Wu & Xue 2002a,b). In these models the 
self-similar entropy distribution of unmodified in- 
tracluster gas is simply truncated below Kc, and 
the gas is allowed to relax back into hydrostatic 
equilibrium. These truncated models agree bet- 
ter with the hotter clusters but not do not pro- 
vide much enhancement over Kc near ~1 keV be- 
cause this simplistic prescription leads to nearly 
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iscntropic cores in the low-mass clusters. Treat- 
ing cooling more realistically produces an entropy 
gradient within 0.1r2oo that declines to zero en- 
tropy at the origin, in the absence of feedback and 
conduction (Voit et al. 2002). 

A more detailed computation of Kc accounting 
for time; despondence of the gas temperature does 
not change this basic result. One can compute the 
total entropy loss Ky^d owing to radiative cooling 
assuming that the gas temperature at a given time 
equals the temperature of the most massive pro- 
genitor halo at that time (see Voit et al. 2003). 
The squares and triangles in Figure 1 show that 
if rad differs little from Kc for a typical merger his- 
tory (see also Oh & Benson 2003). 

Taken as whole, Figure 1 strongly suggests 
that radiative cooling regulates the entropy scale 
that defines the cores of clusters. However, the 
large dispersion of the unbinned data points indi- 
cates that cooling acts in concert with other pro- 
cesses. Much of that dispersion may arise from 
the stochastic nature of mergers, which produce 
significant entropy differences in clusters of simi- 
lar temperature (see Figure 5 of Voit et al. 2003). 
Differing feedback histories may also be responsi- 
ble for some of this dispersion. 

3. Outer Entropy cind Smoothing 

Elevated entropy levels in cluster cores were 
once thought to be evidence for non-gravitational 
energy input into the intracluster medium, but the 
close link between core entropy and the cooling 
threshold now renders that interpretation of core 
entropy somewhat questionable. Nevertheless, 
feedback from galaxy formation still seems nec- 
essary to limit the fraction of condensed baryons 
in clusters. In models of hierarchical structure for- 
mation, a large percentage of the plasma that will 
become the intracluster medium has a cooling time 
less than a Hubble time at z ~ 2-3 (e.g., Voit 
et al. 2002). Because only ~5-10% of baryons 
are observed to be in stars (e.g., Balogh et al. 
2001), feedback from supernovae and active galac- 
tic nuclei presumably intervenes to prevent most 
of this intergalactic gas from condensing. Here we 
present evidence suggesting that this early episode 
of cooling and feedback may couple with subse- 
quent merging and accretion to produce an observ- 
able entropy signature in the outskirts of clusters 



owing to the smoothing effect of galactic feedback 
on the intergalactic medium. 

Hierarchical merging in the absence of cooling 
and feedback tends to produce a nearly self-similar 
entropy distribution (Voit et al. 2003). The char- 
acteristic entropy associated with the baryons in 
a halo of mass M200 is -^"200 = T2oone'^^^ , where 
^200 = GM2oaiJ-mp/2r2ao is the characteristic tem- 
perature of the halo and fie ~ 1-6 x 10~^cm~^ 
is 2000^^ times the mean electron density of the 
universe. According to the simulations of Bryan 
& Voit (2001), the maximum value of intracluster 
entropy is typically « 0.8ii'20O without galaxy for- 
mation, and the median entropy value at the scale 
radius rsoo is « 0.4X200- These scalings do not 
appear to depend significantly on halo mass. 

In contrast, observations imply a considerably 
higher entropy value at r^oo with a clear depen- 
dence on halo mass (Finoguenov et al. 2002; Pon- 
man et al. 2003). Figure 2 shows mean values of 
K{r5oo) inferred by Ponman et al. (2003) from 
the same data set as in Figure 1. Notice that all 
of the binned data points lie above the dotted line 
showing the value of K{r 500) predicted by self- 
similar models calibrated with the simulations of 
Bryan & Voit (2001), indicating that galaxy for- 
mation affects the entire entropy profile of a clus- 
ter, not just the core. The observed entropy boost 
is quite large. A solid line showing the value of 
^200 inferred from Ti^m and the observed M200- 
Tium relation of Sanderson et al. (2003) indicates 
the maximum entropy expected from pure hierar- 
chical accretion. All the data points below ~ 5 keV 
lie above this line by hundreds of keV cm^, rep- 
resenting an entropy excess that is very difficult 
to produce with non-gravitational energy injection 
alone (Ponman et al. 2003). 

Voit et al. (2003) and Ponman et al. (2003) 
independently suggested that this excess entropy 
arises because feedback has smoothed the baryonic 
matter that accretes onto a cluster, thereby reduc- 
ing the characteristic density of gas that passes 
through the accretion shock. Smoothing boosts in- 
tracluster entropy because the entropy of strongly 
shocked gas is oc v^p^'^^^, where Vg is the shock ve- 
locity and pi is the preshock baryon density. The 
scale of Vs is determined by dark-matter dynamics 
that are virtually unaffected by cooling and feed- 
back. Thus, feedback raises the mass-weighted av- 
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Fig. 2. — Entropy at rr,m) as a function of luminosity- 
weighted cluster temperature. Each point with error bars 
shows the mean value of K{r5oa) implied by the density 
and temperature profiles of eight clusters within that tem- 
perature bin (Ponman et al. 2003). The dotted line shows 
the predicted entropy at rsoo for self-similar clusters cali- 
brated with simulation L50-f of Bryan & Voit (2001). The 
solid line shows the entropy scale K200 computed from Ti^^ 
using the M2oo-Tium relation of Sanderson et al. (2003), 
which is the maximum amount of intracluster entropy ex- 
pected from hierarchical axicrection in the absence of non- 
gravitational processes (see Voit et al. 2003). Dashed lines 
connect the maximum entropy values produced by smooth 
accretion onto individual halos, according the accretion his- 
tories of Voit et al. (2003, triangles) and Tozzi & Norman 
(2001, squares). The data suggest a transition from hier- 
archical accretion in the most massive clusters to smooth 
accretion on group scales. 



erage value of v^p^"^^^ because it reduces the mass- 
weighted average vahic of pi (Voit ct al. 2003). 
Notice that a modest amount of preshock entropy 
can be strongly amplified at the accretion shock 
of a deep potential well if that preshock entropy 
is sufficient to drive a large proportion of baryons 
out of the shallower potential wells of the accreting 
dark matter. Ponman ct al. (2003) envisioned this 
smoothing as a thickening of the filaments through 
which matter accretes. Voit et al. (2003) mod- 
eled this entropy increase by developing analyti- 
cal models for entropy production during smooth 
accretion. 

The models of Voit et al. (2003) show that 
smooth, spherically symmetric accretion produces 
an entropy profile whose shape is similar to that 



found in simulations without galaxy formation but 
whose normalization is ~2 times higher in massive 
clusters and up to ~4 times higher in groups. The 
magnitude of this entropy excess depends on the 
mass-accretion rate — objects that accrete most of 
their mass late in time have a higher baryon den- 
sity entering the accretion shock and therefore less 
of an entropy excess. The dashed lines in Figure 2 
indicate the maximum amount of entropy pro- 
duced by smooth accretion, assuming the accre- 
tion histories of Tozzi & Norman (2001; squares) 
and Voit et al. (2003; triangles) for halo masses of 
W^h-^MQ, lOi* /i-^M©, and W^h'^MQ. 

This initial comparison between analytical 
models for smooth accretion and observed entropy 
levels at rsoo suggest that smoothed accretion is 
a plausible explanation for these entropy excesses 
because the data points in Figure 2 approach but 
do not exceed the dashed lines. Furthermore, 
the trend of the data points with Tium suggests 
that accretion of baryons onto groups is smoother 
than accretion of baryons onto clusters. This ap- 
parent transition from nearly smooth accretion 
at < IkeV to purely hierarchical accretion at 
> 10 keV may also explain why the entropy pro- 
files of groups are similar in shape to those of clus- 
ters but have shifted normalizations (e.g., Pratt 
& Arnaud 2003). Curiously, scaling the profiles 
by the same T^/^ factor fotmd in the cores brings 
them into agreement, suggesting that the cooling 
threshold plays a role in setting the smoothing 
scale. However, it is not immediately clear how 
cooling and feedback would couple with accretion 
to produce this scaling. 

4. Discussion 

The entropy of baryons in both the cores and 
outskirts of clusters disagrees with simulations 
that do not include galaxy formation. Thus, the 
entropy profiles of these objects contain valuable 
information about the impact of galaxy formation 
on the intcrgalactic medium and the role of feed- 
back in regulating global star formation. The scal- 
ing of core entropy {Ko_i oc Tj^^^) suggests that it 
is governed by radiative cooling (Ponman et al. 
2003). However, feedback is likely to limit the 
fraction of baryons that eventually cool and form 
stars. 

If feedback is strong enough to eject a sub- 
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stantial percentage of baryons from low-mass ha- 
les, then this smoothing of the baryon distribu- 
tion will affect entropy production during subse- 
quent accretion and mergers. Comparing smooth- 
accretion models with the available data (Fig- 
ure 2) suggests that a large percentage of the mat- 
ter accreting onto groups and low-mass clusters 
must be smoothed in order to explain the observed 
entropy excess at rsoo- If this explanation is cor- 
rect, then smoothing may also have an observable 
impact on the X-ray substructure in and around 
low-mass clusters and the associated filaments. 
Measuring how feedback affects the level of X-ray 
substructure in numerically simulated groups and 
clusters would help observers to quantify the level 
of intergalactic smoothing. 

Only in the highest-mass clusters do entropy 
profiles approach the self-similar form predicted by 
simulations without cooling and feedback. How- 
ever, these clusters presumably formed through 
mergers of lower-mass objects in which cooling 
and feedback had already broken self-similarity. 
We speculate that hierarchical merging restores 
similarity because the self-similar entropy profile 
behaves as a dynamical attractor: As structure 
grows, the deepening potential wells are increas- 
ingly able to overcome baryon smoothing pro- 
duced by early feedback. Accretion onto mas- 
sive clusters is therefore lumpier than accretion 
onto groups, and the resulting entropy profile ap- 
proaches self-similarity when the accreting lumps 
of dark matter are able to retain the majority of 
their associated baryons. 

Many questions about the link between galaxy 
formation and the intraclnster medium remain 
unanswered. Among the most crucial is the 
amount of energy injection implied by the en- 
tropy profiles of clusters. Numerical simulations 
that track entropy production in low-mass halos 
would be very helpful because smoothing allows a 
small amount of entropy injection to be amplified 
with each merger event during the formation of a 
group or cluster. High resolution simulations will 
therefore be needed to assess how early feedback 
propagates through the entire merger hierarchy 
and to characterize how the hottest clusters even- 
tually approach self-similarity. 

The authors thank Greg Bryan for sharing 
his insights into simulated clusters and Alistair 



Sanderson for his contributions to the analysis of 
the observational data. GMV received partial sup- 
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